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Abstract 

Background: Tlie late bliglit patliogen Phytophthora infestans can attack botli potato foliage and tubers. Although 
interaction transcriptome dynamics between potato foliage and various pathogens have been reported, no 
transcriptome study has focused specifically upon how potato tubers respond to pathogen infection. When 
inoculated with P. infestans, tubers of nontransformed 'Russet Burbank' (WT) potato develop late blight disease 
while those of transgenic 'Russet Burbank' line SP221 1 which expresses the potato late blight resistance gene 
RB (Rpi-blbl), do not. We compared transcriptome responses to P. infestans inoculation in tubers of these two lines. 

Results: We demonstrated the practicality of RNA-seq to study tetraploid potato and present the first RNA-seq 
study of potato tuber diseases. A total of 483 million paired end lllumina RNA-seq reads were generated, 
representing the transcription of around 30,000 potato genes. Differentially expressed genes, gene groups and 
ontology bins that exhibited differences between the WT and ^RB lines were identified. P. infestans transcripts, 
including those of known effectors, were also identified. 

Conclusion: Faster and stronger activation of defense related genes, gene groups and ontology bins correlate with 
successful tuber resistance against P. infestans. Our results suggest that the hypersensitive response is likely a 
general form of resistance against the hemibiotrophic P. infestans — even in potato tubers, organs that develop 
below ground. 

Keywords: Solanum tuberosum, Next-generation sequencing (NGS), RNA-seq, Pliytoplitliora infestans, FPKM 
(fragments per kilobase of exon per million mapped reads), JA (jasmonic acid), SA (salicylic acid), ET (Ethylene), 
qRT-PCR 



Background 

Cultivated potato is the worlds third most important 
human food crop and the number one non-grain food 
commodity (FAOSTAT 2010). Unfortunately, potato is 
also host of a broad range of pathogens [1]. Late blight 
disease is caused by Phytophthora infestans^ resulted in 
the Irish Potato Famine in the 1840s, and still today 
results in multi-billion dollar losses worldwide annually. 
P, infestans is a notorious plant destroyer with the 
capacity to attack both potato foliage and tubers. Foliage 
resistance against late blight does not guarantee tuber 
resistance— contrasting disease resistance phenotypes 
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can be evident in comparing foliage and tubers from a 
single genotype [2]. 

Gene RB [Rpi-blbl) [3,4], a disease resistance gene 
cloned from a wild potato species, confers broad- 
spectrum foliar resistance against all major late blight 
pathogen isolates. The gene has been introduced into 
several potato cultivars using transgenic approaches [5]. 
Transgenic lines have been or are currently being tested 
for eventual commercial release in Europe, India, 
Bangladesh, US and other places. Halterman et al. [6] 
reported that the four foliar late blight resistant ^-RB 
transgenic potato lines examined in their study lacked 
statistically significant tuber blight resistance. However, 
in a more extensive survey of +RB transgenic potato 
lines, we identified two transgenic lines with unusually 
high RB transcript levels that were late blight resistant in 
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both foliage and tuber [6]. Thus, RB has the potential 
to function in the tuber and the i?5-potato tuber-P. 
infestans interaction provides a tractable system to study 
how potato tubers defend against plant pathogens. 

Next-generation sequencing (NGS) technologies are 
fast evolving and are transforming biology research [7,8]. 
Genome sequences of potato and P, infestans have been 
published [9,10], making sequencing-based transcrip- 
tome studies (RNA-seq) more accessible to potato late 
blight researchers. RNA-seq is a relatively new approach 
towards study of the transcriptome [11,12]. To our 
knowledge, there has been no RNA-seq study focused 
on potato-microbe interactions. Indeed, one NGS -based 
transcriptome study focused on the interaction between 
potato and P, infestans has been published [13]. How- 
ever, that study utilized the DeepSAGE method, not 
RNA-seq. DeepSAGE differs substantially from RNA-seq 
and researchers in that study relied heavily on assembled 
tags; the newly available genome sequence data were 
mostly not utilized. Furthermore, that and all other pub- 
lished studies focused on potato foliage transcriptome 
dynamics [13-15]; no previous study has reported tran- 
scriptome dynamics of potato tubers in response to 
pathogen attack. 

In this study, we employed RNA-seq to study the 
transcriptome dynamics of potato tuber- P, infestans 
interactions in compatible and incompatible potato 
genotypes. We employed genome-wide sequence data 
from both potato and P, infestans in our analyses. 
Differentially expressed genes and ontology bins were 
identified that distinguish compatible and incompat- 
ible interactions. Transcripts of P, infestans^ including 
those from candidate effectors, were also identified. 
Our study has important implications for potato R 
gene deployment and contributes to scientific under- 
standing of organ-specific defense regulation in plants. 

Results 

The +RB transgenic line exhibits enhanced tuber late 
blight resistance 

In previous field-based evaluations, the ^-RB line exam- 
ined in this study, SP2211, was ranked as the third 
most foliar late blight resistant of 57 transgenic lines 
tested and was rated as "Resistant" to foliar late blight 
[5]. In contrast, the WT line (nontransformed 'Russet 
Burbank') was rated as "Susceptible" to foliar late 
blight [5]. Here, in replicated whole tuber assays 
performed six weeks after harvest, the ^-RB line 
showed no tuber late blight disease after normalization 
to the water-inoculated controls, whereas the WT line 
showed clear tuber late blight disease development 
(Figure 1). Results from two sets of whole tuber experi- 
ments were in agreement (32 of 33 P, infestans inocu- 
lations of WT tubers resulted in disease, none of the 




Figure 1 Gene RB confers enhanced tuber late blight disease 
resistance. Field grown tubers of nontransformed 'Russet Burbank' 
(\/\^ and SP221 1 {+RB), a transformed 'Russet Burbank' line carrying 
the RB transgene, were mechanically wounded, inoculated with 
Phytophthora infestans, and incubated for 1 1 days under conditions 
that favor disease development. Tubers were peeled to allow 
assessment of diseased tissues. WT tubers consistently show robust 
tuber late blight disease development, revealed as darkened tissue 
radiating from inoculation sites. In contrast, +RB tubers displayed no 
tuber late blight disease. Note that brown spots present on +RB 
tubers are in response to mechanical wounding, not late blight 
disease; similar wound response was observed in water-inoculated 
tubers (not shown). 



60 inoculation sites in ^RB tubers developed late blight 
disease). Thus, while WT is susceptible to P. infestans 
infection in both foliage and tuber, the RB transgene 
renders the ^RB line SP2211 resistant to P. infestans 
infection in both foliage and tuber. RB gene transcrip- 
tion in the tubers of the -\-RB line was unaltered by 
inoculation with P. infestans at 0, 24, and 48 hpi 
(p>0.05). 

RNA-seq reads aligned well with the potato reference 
genome sequence 

A total of 36 RNA samples, collected from three bio- 
logical replicates of P. infestans- or water-inoculated 
tuber tissues of two potato genotypes (WT and ^-RB) at 
three time points post inoculation, were subjected to 
RNA-seq. Approximately 483 million paired end reads 
were generated, yielding an average of 13.4 million 
paired end reads per sample. A total of 436.3 million 
Illumina reads (90.3%) passed quality filtering. The ma- 
jority of reads (380.2 million or 78.7% of all reads) could 
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be mapped uniquely to one location within the doubled 
monoploid (DM) potato reference genome sequence [9] 
(Figure 2). An additional 16.6 million reads (3.5%) were 
mapped to multiple locations within the reference genome 
sequence (Figure 2). 

We detected transcription of 29,319 potato genes 
based on cufflinks FPKM (fragments per kilobase of 
exon per million mapped reads) information and gene 
models reported by the Potato Genome Sequencing 
Consortium (PGSC) [9]. Wang et al. [16] demonstrated 
that a data set of 10 million RNA-seq reads represented 
about 80% of annotated chicken genes, concluding that 
RNA-seq is a viable alternative to the microarray for 
transcriptome study. In this study, we detected more 
than 75% of PGSC potato gene models— a percentage 
similar to that reported by Wang et al. [17]. This scale of 
gene detection is higher than reported for any previous 
potato-P. infestans transcriptome study. Importantly, 
19.8% of RNA-seq reads that passed quality filters were 
mapped to regions outside of the current potato reference 
genome gene models (or predicted exons), suggesting 




Figure 2 A majority of RNA-seq reads map uniquely to the 
potato reference genome sequence. A pie chart summarizing 
results of alignment of lllumina RNA-seq reads from 36 water- or 
P. /nfestc/ns-i noculated MMJ and +RB tuber samples to the potato 
reference genome sequence. The grey portion of the chart 
represents RNA-seq reads that failed quality checks and were filtered 
out of our data set (9.6% of all RNA-seq reads). The remaining 90.4% 
of RNA-seq reads passed quality checks. The blue portion of the 
chart represents reads that mapped uniquely to the potato genome 
sequence (78.7% of all RNA-seq reads). The black portion of the 
chart represents reads that mapped to multiple locations of the 
potato genome sequence (3.5%). The yellow portion of the chart 
represents reads that failed to map to the potato genome 
sequence (8.2%), including <0.01% of RNA-seq reads that mapped 
to P. infestans transcripts. 

v J 



further refinement of gene annotations associated with the 
potato reference genome sequence may be warranted. 
Nonetheless, qRT-PCR results correlated well with RNA- 
seq data with an average Pearson correlation coefficient of 
r=0.86 (Additional file 1). Together, our results suggest that 
RNA-seq and qPCR approaches can be cross-validated, 
confirming that current potato gene models [9] are appro- 
priate for functional genomics studies of potato tubers. 

Overall transcriptome dynamics 

A principal component analysis (PC A) of log2 transformed 
FPKM values for 29,319 genes from the 36 RNA-seq 
samples is shown in Figure 3. Together, PCI and PC2 
explained >40% of the total variance of this dataset. Sam- 
ples collected at 0 hpi were distinct from all other samples, 
with PCI providing clear separation of 0 hpi samples from 
24 hpi and 48 hpi samples. PCI also differentiates water- 
and P. infestanS'inoculated samples at 48 hpi. Overall, the 
distribution of -\-RB samples is similar to that of corre- 
sponding WT samples. These results indicate that time 
and treatments {P. infestans- vs. water-inoculation) play a 
greater role than genotype in defining overall transcrip- 
tome dynamics. This results is expected given that the 
^-RB line was created by transformation of the WT line 
with the RE transgene. 

Differentially expressed (DE) genes 

We also identified 2,531 genes that are differentially 
expressed (DE) during at least one of the between time 
transitions in P, infestans 'inoculated samples (Additional 
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Figure 3 Time and treatment have larger influences than 
genotype on overall transcriptome differences. FPKM (Fragment 
per kilobase of exon per million mapped reads) values for 
transcriptome sets from 36 potato tuber samples were subjected to 
PCA using the R statistical package [53]. Blue circles represent P. 
/nfestans-inoculated tuber samples collected at 0 hpi; green circles 
represent water-inoculated samples at 0 hpi; orange circles represent 
P. /nfestans-inoculated samples collected at 24 hpi; purple circles 
represent water-inoculated samples collected at 24 hpi; red circles 
represent P. /nfestans-i noculated samples collected at 48 hpi; black 
circles represent water-inoculated samples collected at 48 hpi. 
Circles containing dots were collected from the transgenic line 
SP221 1 i+RB); circles without dots were collected from 
nontransformed 'Russet Burbank' {WT). Note that tuber samples 
collected at a similar time tend to cluster, regardless of genotype 
of origin. 
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file 2). Transcriptome dynamics of both genotypes in 
P. infestans- compared to water-inoculated samples 
during transitions from 0 hpi to 24 hpi and from 24 hpi 
to 48 hpi were visualized using MA plots (hybridization 
intensity plotted against the fold change in expression; 
Figure 4). Interestingly, at later stages of infection (24 
hpi - 48 hpi), WT tubers exhibited mostly up-regulation 
of genes, whereas the ^-RB line displayed approximately 
equal numbers of up- and down-regulated genes (Figure 4), 
potentially reflecting contrasting rates of pathogen prolifer- 
ation (as indicated by total pathogen reads counts, see 
results sections below) and corresponding host response 
patterns. 

Transcriptome dynamics of between line comparisons 
were also visualized using MA plots (Additional file 3). 
Consistent with PC A which shows that genotype has a 



lesser role (compared to time or treatment) in defining 
overaU transcriptome dynamics, fewer genes were deter- 
mined to be DE during between line comparisons, des- 
pite a trend of increasing numbers of DE genes over 
time (Additional file 3). This suggests that both the WT 
and the ^-RB line use that same set of genes to respond 
to P, infestans attack but that there are significant differ- 
ences in the temporal regulation of these genes. 

Next we employed hierarchical clustering to visualize 
the DE genes identified through between time point 
comparisons. Despite a preponderance of shared DE 
gene regulation patterns, groups of genes that distin- 
guish compatible and incompatible interactions were 
stiU identified (Additional files 3 and 4). These include 
WRKY (PGSC0003DMG400008188) and ethylene re- 
sponse (PGSC0003DMG400025989, ERFl) transcription 




C: +/?BOto 24 hpi D: +/?B 24 to 48 hpi 




-5 0 5 10 -5 0 5 10 

Mean gene expression Mean gene expression 

Figure 4 Between time point transcriptome dynamics reveal different patterns of gene regulation in compatible and incompatible 
potato tuber - Phytophthora infestans interactions. MA plots (hybridization intensity plotted against the fold change in expression) displaying 
between time point transcriptome dynamics in nontransformed 'Russet Burbank' {WT) and transgenic SP221 1 {+RB) lines in response to P. 
infestans. X-axes indicate mean gene expression levels [0.5''(log2(FPKMl)+log2(FPKM2))] across the two selected time points for each comparison. 
Y-axes indicate fold change values [log2(FPKMl/FPKM2), where FPKMl represents the later time point and FPKM2 represents the earlier time 
point] across the selected time points. A: changes in WT transcriptome dynamics from 0 to 24 hpi; B: changes in WT transcriptome dynamics 
from 24 hpi to 48 hpi; C: changes in +RB transcriptome dynamics from 0 hpi to 24 hpi; D: changes in +RB transcriptome dynamics from 24 hpi to 
48 hpi. Within each panel, colored dots represent genes that are significantly (FDR<0.001) differentially regulated among comparisons. Black dots 
represent genes that are not significantly differentially regulated. Note that the WT line shows mostly up-regulation of genes during later stages 
of infection (panel B) while the +RB line displays approximately equal up- and down-regulation of genes (panel D). 
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factors. The ERFl transcription factor gene is a likely 
key element in the integration of J A (jasmonic acid) and 
ET (ethylene) signals for the regulation of defense re- 
sponses [17]. 

The +RB transgenic line has faster and stronger induction 
of DE genes 

A total of 1,102 DE genes were identified when water- 
and P, infestanS'inoculdited samples were compared 
(between treatment comparisons) at the same time 
points (24 hpi or 48 hpi). Of these, 959 were DE only in 
the ^RB line, 34 were DE only in the WT line, and 109 
were DE in both lines. We identified two representative 
hierarchical clusters, comprising a total of 31 genes. 



showing interesting regulation patterns at 48 hpi 
(Figure 5). In WT, these 31 genes were suppressed 
(compared to water-inoculated); in -\-RB, these 31 genes 
were induced (compared to water-inoculated). These 
genes include ethylene response factors {ERF2), cellulose 
synthase, chitinase, and elicitor inducible cytochrome P450. 
Two genes specifically associated with the hypersensitive 
response (HR), PGSC0003DMG400025335 [Hypersensitive- 
induced reaction protein {HIR); Figure 5] and PGSCOO 
03DMG400027473 {Hypersensitive response assisting 
protein) were significantly induced in the ^RB line but 
not in the WT line at 48 hpi. It is important to high- 
light that the ^RB line displayed 903 DE genes at 48 
hpi, accounting for about 82% of the total DE genes 




PGSC0003DMG400002099 Mitochondrial carrier protein 
PGSC0003DMG400001804 Glutathione s-transf erase 
PGSC0003DMG400011015 Receptor protein kinase 
PGSC0003DMG400005339 Conserved gene of unknown function 
PGSC0003DMG400015355 Monooxygenase 
PGSC0003DMG401004358 Dienelactone hydrolase family protein 
PGSC0003DMG400008389 Short chain alcohol dehydrogenase 
PGSC0003DMG400024765 Epoxide hydrolase 
PGSC0003DMG400025335 Hypersensitive-induced reaction protein 
PGSC0003DMG400013398 Cellulose synthase 
PGSC0003DMG400007074 CYP72A58 

PGSC0003DMG400031131 Conserved gene of unknown function 
PGSC0003DMG400026261ATERF-2/ATERF2/ERF2 
PGSC0003DMG400015678 10-hydroxygeraniol oxidoreductase 



PGSC0003DMG400012653 LOB domain-containing protein 41 
PGSC0003DMG400003042 Osmotin 
PGSC0003DMG401000287 Myo-inositol oxygenase 
PGSC0003DMG400001372 B12D 

PGSC0003DMG400031104 Conserved gene of unknown function 
PGSC0003DMG400011842 Chitinase 

PGSC0003DMG400026967 Kunitz-type protease inhibitor KPI-D12 
PGSC0003DMG403020240 Glycerophosphodiester phosphodiesterase 
PGSC0003DMG400015219 Miraculin 
PGSC0003DMG400021966 Laccase 
PGSC0003DMG400026222 Major pollen allergen Dry s 1 
PGSC0003DMG400029745 Actin-101 

PGSC0003DMG400018614 Conserved gene of unknown function 
PGSC0003DMG400013684 Elicitor-inducible cytochrome P450 
PGSC0003DMG400046998S-adenosylmethionine-dependentmethyltransferase 
PGSC0003DMG400020264 Conserved gene of unknown function 
PGSC0003DMG400029736 Gene of unknown function 



Log2fold scales 



Figure 5 Hierarchical clustering of differentially expressed (DE) genes in potato tubers following inoculation with Phytophthora 
infestans. Tubers of nontransformed 'Russet Burbank' {WT) and transgenic SP221 1 {+RB) were inoculated witli P. infestans or water. Tuber samples 
collected 0, 24, and 48 hpi were subjected to RNA-seq, revealing a total of 1,102 DE genes between water- and P. /nfestans-inoculated 
comparisons within the same genotype and the same time. Log2(FPKM_p.inf/FPKM_mock) values were used to cluster 1,102 DE genes 
(FDR<0.001) in Cluster 3.0 [51] using uncentered correlation and the complete linkage method. Results were visualized using Treeview [51]. 
(A) Global visualization of the 1,102 DE genes; (B) A small gene cluster differentially regulated in +RB and WT at 24 and 48 hpi; (C) A small gene 
cluster differentially regulated in +RB and WT only at 48 hpi. Red indicates genes that are up-regulated, green indicates genes that 
are down-regulated. 
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identified during between treatment comparisons 
(Figure 6). These results suggest that the ^-RB Une 
responds to P, infestans attack by rapid differential 
regulation of large sets of genes, whereas the WT line 
is slower to respond (e.g., 48 hpi). 

Ontology bins and gene groups distinguishing 
compatible and incompatible interactions 

Mapman ontology enrichment analysis of DE genes 
(Figure 6) revealed that WRKY transcription factors are 




798 



+RB 48 hpi 



Photosynthesis 
Cell wall 

Lipid metabolism 
Secondary metabolism 
Hormone (BA) 
Hormone (ET) 
Hormone (GA) 
Stress 
Peroxidases 
Transcription.WRKY 
Protein. degradation 
Signaling. receptor.kinases 
Transport 

Figure 6 Tubers of the +RB line have a higher frequency of 
differentially expressed (DE) genes 48 hours post inoculation 
(hpi) with Phytophthora infestans, compared to WT (24 and 48 
hpi) and +RB at 24 hpi. Tuber samples collected 0, 24, and 48 hpi 

were subjected to RNA-seq, revealing a total of 1,102 DE genes 
between water- and P. /nfestans-inoculated comparisons within the 
same genotype and the same time. (A) All 1,102 DE genes were 
analyzed using the "Venn count" function in the limma packages of 
R [53] and results were summarized as a Venn diagram. Red: WT 24 
hpi; yellow: +RB 24 hpi; blue: WT 48 hpi; green: +RB 48 hpi. The 
results show that the +RB line is the main contributor of DE genes 
during water- vs. P. /nfestans-inoculated comparisons. (B) All 1,102 
DE genes were also assigned to a MapMan ontology based on the 
Mercator mapping file (see methods), and subjected to Fisher's exact 
test. Bins in red were significantly up-regulated; bins in blue were 
significantly down-regulated; transcription of bins in white did not 
change significantly. The results indicate that ontology bins 
encompassing ET metabolism and signaling are enriched for DE 
genes in +RB but not in WT at 48 hpi. 



up-regulated in -\-RB but not in WT at 48 hpi. Strikingly, 
13 of the 14 DE (during water- vs. P. infestansAnocu- 
lated comparisons) WRKY transcription factor genes 
show significant induction in the -\-RB line at 48 hpi 
(Additional file 5). In contrast, none of the 14 WRKY 
transcription factor genes showed statistically significant 
(FDR<0.001) induction in the WT line at 48 hpi. Homo- 
logs of several of these WRKY transcription factors have 
been previously reported to be involved in defense 
against Phytophthora spp. and other pathogens [18,19]. 
We hypothesize that DE WRKY transcription factors 
identified in this study are downstream regulatory com- 
ponents of i?5-mediated defenses in potato tubers. 

Very importantly, we discovered that when WT and 
+RB were compared directly, the ^-RB line displays 
stronger transcription of genes within receptor kinase 
and defense related bins, even at 0 hpi (Figure 7 and 
Additional file 6). It is especially relevant to note that 



.^^ .qf .tx*^ 

A*^' A*^' ^' 
<^ <£> <^ 

x^ 

Photosynthesis 
Cell Wall 
Lipid metabolism 
Secondary metabolism 
Hormone (Ethylene ET) 
Hormone (JA) 
Stress. biotic 
Transcription.WRKY 
Protein.targeting 
Protein. degradation. E3 
Signaling. receptor kinases 
Transport 

Figure 7 Stronger activation of defense related genes or gene 
groups correlates with successful tuber resistance against P. 
infestans. Tubers of 'Russet Burbank' {WT) and SP221 1 {+RB) were 
inoculated with P. infestans and water. We compared tine RNA-seq 
FPKM counts for WT and +RB using all 39,031 gene models included 
in the Potato Genome Sequencing Consortium (PGSC) v3 dataset 
(i.e., all genes were included, regardless of whether or not a given 
gene was DE). Genes were grouped into ontology bins using a 
Mapman mapping file. Each column represents a comparison 
between the two genotypes at a defined time point post 
inoculation, as indicated. Bins in blue are transcribed at higher levels 
in WT than in +RB; bins in red are transcribed at higher levels in +RB 
than in WT; bins in white did not significantly differ in transcript 
levels between WT and +RB. Results indicate that faster and stronger 
activation of defense bins, most notably biotic stress response and 
receptor kinase bins, occurred in tubers of the tuber late blight 
resistant +RB line. 
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many of the bins commonly associated with plant 
defense, including receptor kinases and PR (pathogen- 
esis-related) protein genes are more highly expressed 
in the ^-RB line. Interestingly, WT actually has much 
higher transcription of similar defense-related bins at 
24 hpi (but not 48 hpi) if an indirect comparison 
(adjusting FPKM values of P, infestans -moculsited sam- 
ples by FPKM values of water-inoculated samples of 
the same genotype) was made, indicating that WT and 
-\-RB probably share large sets of defense components. 
Thus, WT responds to P. infestans by up-regulation of 
sets of defense response genes that are constitutively 
transcribed at higher levels in the ^-RB line even in the 
absence of the pathogen (i.e., 0 hpi). But slower activa- 
tion of defense related genes in WT is insufficient to 
prevent disease development, resulting in a tuber 
blight susceptible phenotype. These results were cor- 
roborated by an independent pilot RNA-seq study 
(lUumina GAIIx single end reads) using field grown 
tubers of the same genotypes from a different year 
(Additional file 7). Collectively, these findings suggest 
that faster and stronger expression of defense related 
genes plays a role in enhanced disease resistance in 
the ^-RB line. 



RNA-seq reads mapping to Phytophthora infestans 
reference transcripts 

A total of 17,353 paired end sequence reads were 
mapped to the P. infestans reference transcript set [10], 
representing the transcription of over 4,600 genes. 
Ninety-six P, infestans transcripts are predicted to 
encode RxLR effectors [10]. 

Throughout the infection time course, a trend of 
increasing mappable P, infestans reads is evident for 
WT samples. This increase of mappable reads was not 
seen in the ^-RB line from 24 hpi to 48 hpi (Additional 
file 8) and a majority (>82%) of all mapped reads 
originated from infected WT potato tuber samples 
collected at 48 hpi. This is in agreement with tuber late 
blight disease development observed in WT but not 
-\-RB tubers. In total, a majority of all P, infestans reads 
(>95%) were from infected WT potato tuber samples; 
less than 5% of the reads were from infected ^-RB 
potato tuber samples (Additional file 9). 

Class I and II ipiO proteins, members of the RxLR family 
of putative effectors, induce avirulence responses in potato 
lines containing RB [20-22]. After inoculation (24, 48 hpi), 
class I ipiO transcripts were detected (though in low 
abundance) in five of six tuber samples collected from the 
P, infestansAnoculated WT line but from no samples 
derived from the -\-RB line (Additional file 10). This is 
probably due to suppression of P. infestans proliferation in 
the -\-RB line. 



Discussion 

Faster and stronger activation of defense related genes, 
gene groups and ontology bins correlates with successful 
tuber defense 

Direct comparison between WT and -\-RB lines at 0 hpi 
revealed a total of 11 DE genes (Additional file 3). 
Seven of the eleven DE genes were more highly tran- 
scribed in ^RB, These include Carbonic anhydrase 
(PGSC0003DMG400006956), Aquaglyceroporin (PCS 
C0003DMG400009604), Cytochrome P450 (PGSC0003 
DMG400015185), Malic enzyme (PGSC0003DMG401 
026923), and 1-aminocyclopropane-l-carboxylate oxidase 
(ACO) (PGSC0003DMG401026923). Interestingly, these 
genes or homologs of these genes were previously reported 
to be involved in basal or pathogen-induced defense re- 
sponses in plants [14,23,24]. The four remaining DE genes 
[protoporphyrinogen oxidase (PGSC0003DMG400001117), 
phenylcoumaran benzylic ether reductase (PGSC0003D 
MG400003691), DNA binding protein (PGSC0003DMG 
400003902), Major Latex (PGSC0003DMG400008811)] are 
transcribed at higher levels in WT than in the -\-RB line. 
Each of these genes has also been associated with plant 
defense [25-28]. Together, our results demonstrate that the 
-\-RB line differs from WT in the constitutive transcriptional 
regulation of defense related genes. 

GO enrichment analysis of DE genes identified during 
between treatment comparisons (Figure 6) suggests that 
+RB has faster response to pathogen attack, as evidenced 
by higher and statistically significant induction of ontology 
bins commonly associated with plant defense (e.g., ET, 
WRKY, Signaling receptor kinases) at 48 hpi. One notable 
example is WRKY transcription factors (Additional file 5), 
with 13 out of 14 WRKY transcription factors significantly 
induced in +RB. These same transcription factors are also 
induced, but at nonsignificant levels, in the WT line, 
documenting that more potent and rapid activation of 
defense related genes, rather than novel resistance re- 
sponses per se, might play a role in enhanced tuber resist- 
ance in the -\-RB line. 

Furthermore, when a more holistic approach (using all 
available data points: 39,031 gene transcript values per 
sample) was employed to identify differentially regulated 
ontology bins, various defense related ontology bins 
were discovered (Figure 7). ^-RB potato tubers displayed 
faster and stronger transcription of genes within bins com- 
monly associated with plant defense such as receptor 
kinases and pathogenesis -related (PJ^ genes (Figure 7). To- 
gether, these results suggest that the ^-RB line is tuber late 
blight resistant due to faster and stronger activation of 
defense related components. 

Cao et al. [29] found that increased transcription of 
the rice Xa3 resistance gene correlated with enhanced 
expression of defense-responsive genes and an enlarged 
resistance spectrum to Xanthomonas oryzae pv. oryzae 
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(Xoo). Previously we demonstrated that higher transgene 
RB copy numbers correlated with both higher RB 
transcript levels and enhanced late blight resistance in 
the foliage [5]. The ^-RB line employed in this study, 
SP2211, displayed the highest tuber transgene transcrip- 
tion levels among 11 transgenic lines examined [6]. 
Thus, the unusually high levels of RB gene transcription 
in the tubers of SP2211 correlates with faster and stron- 
ger activation of defense related components (Figures 6, 7). 
The faster and stronger up-regulation of these defense- 
related components very likely contributed to the observed 
successful defense against P, infestans (Figure 1) in this 
line. One interesting observation is that WT actually has 
much higher transcription of similar defense-related bins 
at 24 hpi (but not 48 hpi) in indirect (adjusting FPKM 
values of P, mfestansAnoculated samples by FPKM values 
of water-inoculated samples of the same genotype) 
comparisons. Collectively, our data support a model of 
tuber blight resistance resulting as a function of faster and 
stronger expression of defense related genes due to high 
RB transcript levels. 

HR is likely a general form of resistance to P. infestans, 
even in potato tubers 

Prevailing scientific thought suggests that the phyllosphere 
frequently encounters biotrophic or hemibiotrophic patho- 
gens, against which the HR is effective. In contrast, roots 
are more likely to encounter necrotrophic microbes that 
would theoretically benefit from the HR [30]. Consistently, 
working in Arabidopsis, Hermanns et al. [31] showed that 
HR is present in above ground incompatible host-pathogen 
interactions, but absent during the same host-pathogen 
interaction in the roots— despite R gene transcription in 
both leaves and roots. These authors concluded that R 
gene function is modulated in an organ-specific but mech- 
anistically poorly defined manner for the purpose of sup- 
pressing HR in the roots. Observations of differential late 
blight resistance levels in potato foliage and tubers within a 
single genotype [2] suggest that R gene function in potato 
may also be modulated in an organ-specific manner. But 
the potato tuber, although it develops below ground, is a 
modified stem, not a root. Previous studies [32,33] have 
suggested that HR plays a role in tuber defense. The identi- 
fication of SP2211, a -\-RB transgenic line with both foliar 
and tuber late blight resistance provided opportunity to 
examine the role of HR in defending non-root organs 
against pathogen attack. 

P. infestans is a hemibiotrophic pathogen and Kamoun 
et al. [34] argued that HR is likely a general form of re- 
sistance against Phytophthora pathogens. In agreement, 
Chen and Halterman [35] demonstrated that RB triggers 
an HR in potato foliage challenged with P, infestans. In 
examining i?5-mediated tuber responses to P, infestans^ 
we identified DE genes associated with the HR (potato 



HIR and hypersensitive response assisting genes) that 
are highly up-regulated in pathogen-inoculated tuber 
samples of SP2211 compared to water-inoculated tuber 
samples at 48 hpi. Previous studies have shown that 
the pepper homolog of HIR is a positive regulator of 
hypersensitive cell death [36,37]. Qi et al. [38] also dem- 
onstrated that the NBS-LRR protein RPS2 forms com- 
plexes with AtHIR proteins in Arabidopsis and tobacco 
plants, providing mechanistic insight into R protein 
function and the HR. We hypothesize that RB and HIR 
might form similar complexes in potato, a topic that 
warrants further experimentation. Rapid production of 
reactive oxygen species by plants is a hallmark of patho- 
gen recognition and correlates with the HR. Working 
with disks cut from potato tubers that carry the Rl late 
blight resistance gene, Doke [32] demonstrated that in- 
compatible but not compatible races of P, infestans trig- 
gered production of reactive oxygen species. In our 
study, a gene (PGSC0003DMG400024754) encoding re- 
spiratory burst oxidase homolog protein B (NADPH oxi- 
dase RBOHB) was induced (based on between treatment 
comparisons) 6.6 fold in the ^-RB line but not in the WT 
line 48 hours after P, infestans -inoculdition. This gene is 
known to be involved in the massive phase II oxidative 
burst induced in potato by pathogen infection [39] 
and offers further support that i?5-mediated resistance 
to tuber late blight likely entails an HR or HR-like 
phenomenon. 

Interestingly, HR is commonly associated with SA 
(salicylic acid) mediated responses. But our ontology bin 
analysis didn't reveal SA metabolism as a differentiating 
factor for the tuber blight resistant and susceptible lines. 
Instead, our study reveals that the ET (ethylene) bin is 
predominantly associated with the incompatible inter- 
action. Leon-Reyes et al. [40] reported that the antagon- 
istic relationship between SA and J A [41] only exists 
when there is no strong production of ET. Nunez- 
Pastrana et al. [42] reported that ethylene but neither SA 
nor Methyl JA (jasmonic acid) induces resistance re- 
sponse against Phytophthora capsici in pepper. Exten- 
sive hormonal crosstalk [43] during plant responses to 
pathogens is a topic of extensive ongoing research and 
future results will likely yield new insights that will 
allow fine-tuning of our understanding of potato tuber- 
pathogen interactions and organ-specific defense responses 
in plants. 

Conclusion 

We presented the first RNA-seq (or transcriptome 
dynamics) study focused on potato tuber responses to 
pathogen attack. We mapped approximately 400 million 
RNA-seq reads onto the recently published potato refer- 
ence genome sequence, documenting the utility of RNA- 
seq for biological study of tetraploid cultivated potato. 
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We identified sets of DE genes that distinguish resistant 
and susceptible Unes. Our data suggest that potent regu- 
lation of defense genes and gene groups or ontology bins 
(e.g., HIR and WRKY transcription factors) plays a role 
in 7?5-mediated tuber defense. In particular, faster and 
stronger expression of defense related genes correlates 
with enhanced tuber blight resistance in the ^-RB trans- 
genic line. In agreement with Kamoun et al. [34] and 
Doke [32], our data suggest that HR is likely a general 
form of resistance against late blight even in the potato 
tuber, an organ that develops below ground. 

Methods 

Plant material and RNA preparation 

The tuber late blight susceptible nontransformed 'Russet 
Burbank' (WT, provided by Dr. Carl Rosen, UMN) and 
SP2211, a tuber late blight resistant transformed 'Russet 
Burbank' line carrying the RB transgene {^-RB [5]) were 
examined in this study. Tubers were produced under 
standard cultural practices at the University of Minne- 
sota Sand Plain Research Farm (Becker, MN). Tubers of 
each line were harvested and held for three days at room 
temperature before storage for six weeks at 11-13 
degrees Celsius. P. infestans US8 isolate US940480 [3] 
was maintained on Rye A medium [44]. Prior to inocula- 
tions, sporangia were harvested from plates by physical 
scraping into distilled water. Inoculum was adjusted to 
6.75*10^ sporangia/ml. Prepared inoculum was incu- 
bated for 1 hour at 4 degrees Celsius and then at room 
temperature for 30 minutes. Six week old WT and ^-RB 
tubers were inoculated using a modified whole tuber 
assay [6], as detailed below. 

For RNA extractions and RNA-Seq, six tubers of each 
genotype were randomly selected, washed with deionized 
water, and allowed to air dry at room temperature 
for 24 hours. Tubers were wounded [0.2 cm*0.3 cm 
(depth*diameter)] at six sites spaced uniformly (~3 cm 
apart) across the tuber surface. Each wound was inocu- 
lated with 10 (il sporangial suspension (3 tubers per 
genotype) or water (mock treatment; 3 tubers per geno- 
type). Inoculated tubers were stored in tightly sealed 
dark boxes under high (-95%) humidity and at room 
temperature for 72 hours (during which time all RNA- 
seq samples were collected as specified below) and then 
moved to 11-13 Celsius for disease development and 
phenotyping at 11 days post inoculation. From each of 
three replicate tubers per genotype, tuber tissue (0.7 cm* 
(0.5-0.8)cm) was collected using a cork borer from an 
inoculation site at time point 0 (pre-inoculation) and at 
6, 12, 24, and 48 hours post inoculation (hpi). The sam- 
pled tissue includes cells from the periderm, cortex and 
medulla layers. The sixth inoculation site on each tuber 
was left intact and the tubers were phenotyped for late 
blight disease development 11 days later by multiplying 



measured length, width, and depth of disease lesions. 
Collected tissue samples (one tuber core per replicate 
tuber for each time point) were immediately frozen in li- 
quid nitrogen and stored at -80 degrees Celsius. To re- 
duce possible experimental variation, tissues collected 
for RNA-seq originated from a common experiment in- 
volving a single P, infestans inoculum preparation or 
water. In this study, each tuber is considered a biological 
replicate of the corresponding genotype by treatment 
combination. In total, 36 tissue samples from the two 
plant genotypes (WT and ^-RB) x three time points (0, 
24, and 48 hpi) x two inocula {P, infestans or water) x 
three replicates were employed for RNA extraction. 
[Note: We collected tissue at 0, 6, 12, 24, and 48 hpi, but 
only 0, 24, and 48 hpi samples were subjected to RNA- 
seq]. Total RNA was extracted from frozen tissue using 
the SV Total RNA Isolation System (Promega Corpor- 
ation, Madison, WI) according to manufacturers in- 
structions. The quantity and quality of RNA samples 
were assessed using a Nanodrop 1000 machine (Thermo 
Fisher Scientific Inc., Wilmington, DE). High quality 
total RNA (5ug, 100 ng/(il) samples were sent to the 
University of Minnesota BioMedical Genomics Center 
(BMGC) for RNA-seq library preparation using the 
TruSeq SBS Kit (50 Cycles) (Illumina Inc., San Diego, 
CA) and paired end sequencing using an Illumina Hi-Seq 
2000 machine (Illumina). 

A separate, additional set of tubers (nine for WT and 
eighteen for ^-RB) were simultaneously subjected to the 
whole tuber assay (using the same inoculum and proto- 
col described above), but without tissue collection, for 
the purposes of disease phenotyping. 

RNA-seq reads mapping and DE genes clustering 

RNA-seq reads were quality filtered using SolexaQA 
packages [45] with default parameters and a length filter 
of greater than 27 bp for both ends of each read 
pair. Sequence data have been submitted to the NCBI 
Sequence Read Archive [accession number SRP022916]. 
Quality filtered RNA-seq reads were analyzed using the 
"Tuxedo Suite" software packages [46]: Bowtie [47] 
vO.12.7, Tophat [48] vl.3.2. Cufflinks [49] vl.1.0 and a 
reference genome of potato [9]. We suppressed Tophat 
[48] from identif)^ing novel junctions. A total of eleven 
pair-wise comparisons [between time (4), between line 
(3), between treatment (4)] were made using Cuffdiff of 
the Cufflinks software packages. Differentially expressed 
(DE) genes are those genes that showed significantly 
(FDR adjusted p-value <0.001) different transcript levels 
among comparisons. False discovery rate correction (FDR) 
was done using the Benjamini Hochberg method [50]. 

Quality filtered RNA-seq reads were also mapped to 
reference transcript sequences of P, infestans [10] and 
different classes (see Champouret et al. [22] for details 
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of classification) of ipiO (potential cognate avr effector 
of RB) gene sequences. The mapping was done using 
Bowtie [47] v0.127. Cufflinks software suites were not 
used in this analysis as underlying normalization assump- 
tions were not met. Potato gene expression fold change 
values were used to cluster and partition genes into groups 
using hierarchical clustering and the complete linkage 
method in Cluster 3.0 [51] and visualized in Treeview 
software [52] using adjusted pixel settings (threshold 1.0). 

The UNIX command line options for Tophat and 
Cuffdiff were as follows: "tophat -G genes.gff -o 
sample.out — no-novel-juncs -r 100 genome.fasta read, 
left.fq read.right.fq; cuffdiff -N -u -O cuffdiff.out 
genes.gff l.bam,2.bam,3.bam 4.bam,5.bam,6.bam". The 
UNIX command line options for Bowtie were as fol- 
lows: "Bowtie -a -v 2 —best -M 1 — fr ref_bowtie_build 
-1 readl.fq -2 read2.fq". Custom perl scripts were 
used to parse the mapping/alignment results into 
tabular formats. R statistical software [53] was used to 
generate various plots. 

Functional assignment and MapMan analysis of potato 
genes 

To assign potato genes into functional categories (bins), 
we adopted MapMan ontology [54]. The Mercator (Max 
Planck Institute, Germany, Dr. Mark Lohse, personal 
communication) annotation pipeline was used to assign 
potato genes into functional bins by searching a variety 
of reference databases. The mapping file was manually 
curated to remove suspicious assignments. Water- and 
P, inf es tans Anoculsited samples were compared using 
Fisher s exact test to identify ontology bins enriched with 
significantly up or down regulated DE (FDR <0.001) genes. 
WT and -\-RB lines were compared using the Wilcoxon 
rank-sum test to identify functional bins that behaved 
differently. Resulting p-values from Fisher s exact tests or 
Wilcoxon rank-sum tests were subjected to multiple test 
correction using the Benjamini Hochberg method [50]. 

Real-time quantitative RT-PCR 

For determining the transcript levels of the transgene 
RB, we used a previously established qRT-PCR protocol 
[5,55]. For RNA-seq validation. Primer Express 3.0 
(Applied Biosystems, Foster City, CA, USA) was used 
with default parameters to generate primer pairs for 
selected transcripts (Additional file 11). Total RNA 
(150 ng, 30 ng/ul) from diluted stocks of the same RNA 
that was subjected to RNA-seq was used in each reverse 
transcription reaction using the Superscript III First 
Strand Synthesis Kit (Life Technologies Inc., Carlsband 
CA) according to manufacturer s instructions. All qPCRs 
were performed using the Power SYBR Green Master 
Mix (Life Technologies) and an ABI 7500 Real Time 
System (Applied Biosystems). Target gene transcript 



levels were normalized to EFla [56] using 2' values. 
Normalized gene expression levels were compared with 
RNA-seq FPKM values derived from Cufflinks. 

Additional files 



Additional file 1: RNA-seq FPKM and qPCR correlation. First column 
is potato tuber sample/IDs. Remair^ing columr^s ir^dicate qPCR results ar^d 
RNA-seq FPKM values for each potato gene, as indicated. The per gene 
and average correlations between RNA-seq and qPCR are listed at the 
bottom of the page. 

Additional file 2: Fold change values for transcript levels of each of 
the 39,301 potato genes in each of the 20 comparisons. 

Additional file 3: MA plots (hybridization intensity plotted against 
the fold change in expression) of between line comparisons. Each 
X-axis indicates mean gene expression levels [log2(FPKMl)+log2(FPKM2)] 
across the two time points selected for comparison. Each y-axis indicates 
fold change values [log2(FPKI\/ll/FPKM2)]. (A) WT compared to +RB at 0 hpi; 
(B): WT compared to +RB at 24 hpi; (C): WT compared to +RB at 48 hpi. 

Additional file 4: Hierarchical clustering and Treeview visualization 
of 1,767 DE genes (determined based on between time point 
comparisons; Note: genes that are DE in between time point 
comparisons in water-inoculated samples were exclude from this 
analysis) show that DE gene regulation patterns in compatible and 
incompatible interactions are predominantly similar. Each column 
represents a comparison between two time points. Column one: 0 hpi to 
24 hpi in V\^; column two: 0 hpi to 24 hpi in +RB; column three: 24 hpi 
to 48 hpi in V\^; column four: 24 hpi to 48 hpi in +RB. Red indicates 
up-regulation, green indicates down regulation. Left panel: the overall 
pattern of the 1,767 genes. Right panel: Magnified images of small gene 
clusters. 

Additional file 5: Transcription of WRKY genes is highly induced in 

+RB but not WT. Column one and two are gene ID and PGSC 
annotation descriptions. Columns three to six are log2 fold change values 
derived from between treatment comparisons {P. infestans- vs. water- 
inoculated). Values highlighted in red are statistically different between 
the +RB and WT lines. 

Additional file 6: Tubers of 'Russet Burbank' (WT) and SP221 1 {+RB) 
were inoculated with P. infestans and water. We compared the RNA- 
seq FPKM counts for V\^ and +RB using all 39,031 gene models included 
in the Potato Genome Sequencing Consortium (PGSC) v3 dataset (i.e., all 
genes were included, regardless of whether or not a given gene was DE). 
Genes were grouped into ontology bins using a Mapman mapping file. 
Each column represents a comparison between the two genotypes at a 
defined time point post inoculation, as indicated. Bins in blue are 
transcribed at higher levels in V\^than in +RB; bins in red are transcribed 
at higher levels in +RB than in WT; bins in white did not significantly 
differ in transcript levels between V\^ and +RB. Results indicate that 
stronger activation of defense bins, including stress responses and 
receptor kinases, occurred in +RB (the tuber blight resistant line). 

Additional file 7: In a pilot study, WT and +RB samples at 0, 24 and 
48 hours post P. infestans inoculation were collected from three 
bio-reps and pooled into a single composite sample for RNA-seq. A 

total of 146.6 million single end lllumina reads (51 bp) were filtered and 
mapped to the reference potato genome using SolexaQA and Tuxedo 
software suite packages [9]. Cuffdiff was used to generate log2 
transformed fold change values for each between genotype comparisons 
at 0, 24 and 48 hpi. Mapman analyses and Wilcoxon rank sum tests were 
performed (see methods). Columns one, two, and three represent +RB vs. 
V\^ comparisons at 0, 24 and 48 hpi, respectively. Blue bins show higher 
transcription in V\^; red bins show higher transcription in +RB. Note that 
+RB has faster and stronger activation of defense related bins (stress and 
receptor kinases) at 24 and 48 hpi. 

Additional file 8: Log2 transformed total read counts that mapped 
to Phytophthora infestans transcripts. The X-axis indicates different 
time points (0, 24, 48 hpi) post P. infestans inoculation. The Y-axis 
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indicates log2 transformed total mapped reads count. Results indicate an 
increase in P. infestans RNA-seq reads in the V\^ but not the +RB line 
over time. 

Additional file 9: Bowtie RNA-seq read mapping results for 4,634 
Phytophthora infestans genes. The first column contains gene ID, the 
second column indicates whether the gene is a predicted RxLR effector 
gene, each remaining column contains mapping summary results for 
each of the 36 potato tuber samples. 

Additional file 10: Bowtie RNA-seq reads mapping results for ipiO 
genes. The first and second columns are ref id and short IDs for the ipiO 
genes adopted from Champouret et al. [22]. Each remaining column 
contains a bowtie mapping summary for each sample following tuber 
inoculation with Phytophthora infestans. 

Additional file 11: qRT-PCR primer sequences for 20 selected 
potato genes. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

LLG was responsible of all bench work and informatics analyses of RNA-seq 
reads filtering, mapping, DE genes detection, and ontology regulation 
analysis. ZJT provided some informatics support, especially in relation to a 
pilot RNA-seq project. BPM identified transgenic line SP221 1 i+RB), pioneered 
tuber resistance phenotyping, and developed molecular assays for transgene 
transcription measurement. JMB provided the conceptual impetus for this 
work, suggestions on all aspects of experimental design and data analysis, 
and secured funding for this project. LLG wrote this manuscript and ZJT, 
BPM, and JMB provided comments and suggestions to improve the 
manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

The work was funded by USDA NIFA, EPA, MN area II Potato Research and 
Promotion Council and a Thesis Research Grant to LLG from the Graduate 
School, University of Minnesota. The authors thank Dr. Mark Lohse of Max 
Planck Institute (Germany) for generating a draft PGSC mapping file for 
MapMan analysis in potato. The authors thank Dr. Fumiaki Katagiri for 
thoughtful comments on statistical analysis and draft versions of our 
manuscript. Computing resources from the Minnesota Supercomputing 
Institute at the University of Minnesota are greatly appreciated. 

Author details 

^Department of Plant Pathology, University of Minnesota, Saint Paul, MN 
55108, USA. ^Biomedical Statistics and Informatics Division, Mayo Clinic, 
Rochester, MN 55905, USA. 

Received: 17 December 2012 Accepted: 17 May 2013 
Published: 23 May 2013 

References 

1. Stevenson WR, Loria R, Frank GD, Weingartner DP: Compedium of potato 
diseases. Saint Paul, MN, USA: APS Press; 2001. 

2. Kirk WW, Felcher KJ, Douches DS, Niemira BA, Hammerschmidt R: Susceptibility 
of potato {Solonum tuberosum L.) foliage and tubers to the US8 genotype of 
Phytophttiora infestans. Am J Potato Res 2001 , 78(4):31 9-322. 

3. Song JQ, Bradeen JM, Naess SK, Raasch JA, Wielgus SM, Haberlach GT Liu J, 
Kuang HH, Austin-Phillips S, Buell CR, et al: Gene RB cloned from Soianum 
bulbocostonum confers broad spectrum resistance to potato late blight. 
Proc Natl Acad Sci USA 2003, 1 00(1 6):91 28-91 33. 

4. van der Vossen E, Sikkema A, Hekkert BTL, Gros J, Stevens P, Muskens M, 
Wouters D, Pereira A, Stiekema W, Allefs S: An ancient R gene from the 
wild potato species Solonum bulbocostonum confers broad-spectrum 
resistance to Phytophthoro infestons in cultivated potato and tomato. 
Plant J 2003, 36(6):867-882. 

5. Bradeen JM, Lorizzo M, Mollov DS, Raasch J, Kramer LC, Millett BP, Austin- 
Phillips S, Jiang JM, Carputo D: Higher copy numbers of the potato RB 
transgene correspond to enhanced transcript and late blight resistance 
levels. Mol Plant Microbe Interact 2009, 22(4):437-446. 



6. Millett BP: Characterization of transcription and function of foliar blight 
resistance gene RB in foliage and tubers of transgenic cultivated potato at 
different physiological ages. University of Minnesota: PhD Thesis; 2008:178. 

7. Morozova 0, Marra MA: Applications of next-generation sequencing 
technologies in functional genomics. Genonnics 2008, 92(5):255-264. 

8. Metzker ML: Applications of next generation sequencing sequencing 
technologies - the next generation. Nat Rev Genet 2010, 11(1):31-46. 

9. Xu X, Pan SK, Cheng SF, Zhang B, Mu DS, Ni PX, Zhang GY, Yang S, Li RQ, 
Wang J, et al: Genome sequence and analysis of the tuber crop potato. 
Nature 201 1, 475(7355):189-U194. 

10. Haas BJ, Kamoun S, Zody MC, Jiang RHY, Handsaker RE, Cano LM, Grabherr 
M, Kodira CD, Raffaele S, Torto-Alalibo T, et al: Genome sequence and 
analysis of the Irish potato famine pathogen Phytophthoro infestons. 
Nature 2009, 461(7262):393-398. 

11. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and 
quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 2008, 
5(7):621-628. 

12. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: An 
assessment of technical reproducibility and comparison with gene 
expression arrays. Genome Res 2008, 18(9):1509-1517. 

13. Gyetvai G, Sonderkaer M, Gobel U, Basekow R, Ballvora A, Imhoff M, Kersten 
B, Nielsen KL, Gebhardt C: The transcriptome of compatible and 
incompatible interactions of potato {Solonum tuberosum) with 
Phytophthoro infestons revealed by DeepSAGE analysis. PLoS One 2012, 
7(2):e31526. 

14. Restrepo S, Myers KL, del Pozo 0, Martin GB, Hart AL, Buell CR, Fry WE, 
Smart CD: Gene profiling of a compatible interaction between 
Phytophthoro infestons and Solonum tuberosum suggests a role for 
carbonic anhydrase. Mol Plant Microbe Interact 2005, 18(9):91 3-922. 

1 5. Birch PRJ, Avrova AO, Duncan JM, Lyon GD, Toth RL: Isolation of potato genes 
that are induced during an early stage of the hypersensitive response to 
Phytophthoro infestons. Mol Plant Microbe Interact 1 999, 1 2(4):356-361 . 

16. Wang Y, Ghaffari N, Johnson CD, Braga-Neto UM, Wang H, Chen R, Zhou HJ: 
Evaluation of the coverage and depth of transcriptome by RNA-Seq in 
chickens. BMC Bioinforma 2011, 12(suppl 10):S5. 

17. Lorenzo 0, Piqueras R, Sanchez-Serrano JJ, Solano R: Ethylene response 
factorl integrates signals from ethylene and jasmonate pathways in 
plant defense. Plant Cell 2003, 15(1):165-178. 

18. Nakazawa-Ueji YE, Nunez-Pastrana R, Souza-Perera RA, Santana-Buzzy N, 
Zuniga-Aguilar JJ: Mycelium homogenates from a virulent strain of 
Phytophthoro copsici promote a defence-related response in cell 
suspensions from Copsicum chinense. Eur J Plant Pathol 2010, 1 26(3):403-41 5. 

19. Katou S, Yoshioka H, Kawakita K, Rowland 0, Jones JDG, Mori H, Doke N: 
Involvement of PPS3 phosphorylated by elicitor-responsive mitogen- 
activated protein kinases in the regulation of plant cell death. 

Plant Physiol 2005, 139(4):1914-1926. 

20. Pieterse CMJ, Vanwest P, Verbakel HM, Brasse P, Vandenbergvelthuis GCM, 
Covers F: Structure and genomic organization of the ipib and ipio gene 
clusters of Phytophthoro-infestons. Gene 1 994, 1 38(1 -2):67-77. 

21. VIeeshouwers VGAA, Rietman H, Krenek P, Champouret N, Young C, Oh S-K, 
Wang M, Bouwmeester K, Vosman B, Visser RGF, et al: Effector genomics 
accelerates discovery and functional profiling of potato disease 
resistance and Phytophthoro infestons avirulence genes. PLoS One 2008, 
3(8):e2875. 

22. Champouret N, Bouwmeester K, Rietman H, van der Lee T, Maliepaard C, 
Heupink A, van de Vondervoort PJI, Jacobsen E, Visser RGF, van der Vossen LAG, 
et al: Phytophthoro infestons isolates lacking class I ipiO variants are virulent 
on Rpi-blbl potato. Mol Plant Microbe Interact 2009, 22(1 2):1 535-1 545. 

23. Yu MM, Shen L, Zhang AJ, Sheng JP: Methyl jasmonate-induced defense 
responses are associated with elevation of 1-aminocyclopropane-1- 
carboxylate oxidase in Lycopersicon esculentum fruit. J Plant Physiol 201 1, 
1 68(1 5):1 820-1 827. 

24. Ishihama N, Yamada R, Yoshioka M, Katou S, Yoshioka H: Phosphorylation 
of the Nicotiono benthomiono WRKY8 transcription factor by MARK 
functions in the defense response. Plant Cell 201 1 , 23(3):1 1 53-1 1 70. 

25. Molina A, Volrath S, Guyer D, Maleck K, Ryals J, Ward E: Inhibition of 
protoporphyrinogen oxidase expression in arabidopsis causes a lesion- 
mimic phenotype that induces systemic acquired resistance. Plant J 1999, 
17(6):667-678. 

26. Gang DR, Kasahara H, Xia ZQ, Vander Mijnsbrugge K, Bauw G, Boerjan W, 
Van Montagu M, Davin LB, Lewis NG: Evolution of plant defense 



Gao et al. BMC Genomics 201 3, 14:340 Page 1 2 of 1 2 

http://www.bionnedcentral.conn/1 471 -21 64/1 4/340 



mechanisms - relationships of phenylcoumaran benzylic ether 
reductases to pinoresinol-lariciresinol and isoflavone reductases. J Biol 
Chem 1999, 274(1 1):75 16-7527. 

27. Koffhage U, Trezzini GF, Meier I, Hahlbrock K, Somssich IE: Plant 
homeodomain protein involved in transcriptional regulation of a 
pathogen defense-related gene. Plant Cell 1 994, 6(5):695-708. 

28. Konno K: Plant latex and other exudates as plant defense systems: roles 
of various defense chemicals and proteins contained therein. 
Phytochemistry 201 1 , 72(1 3):1 5 1 0-1 530. 

29. Cao YL, Ding XH, Cai M, Zhao J, Lin YJ, Li XH, Xu CG, Wang SP: Expression 
pattern of a rice disease resistance gene Xa3/Xa26 is differentially 
regulated by the genetic backgrounds and developmental stages that 
influence its function. Genetics 2007, 177(1):523-533. 

30. Glazebrook J: Contrasting mechanisms of defense against biotrophic and 
necrotrophic pathogens. Annu Rev Phytopothol 2005, 43:205-227. 

31. Hermanns M, Slusarenl<o AJ, Schlaich NL: Organ-specificity in a plant 
disease is determined independently of R gene signaling. Mol Plant 
Microbe Interact 2003, 16(9):752-759. 

32. Doke N: Involvement of superoxide anion generation in the 
hypersensitive response of potato - tuber tissues to infection with an 
incompatible race of Phytophthora infestans and to the hyphal wall 
components. Physiological Plant Pathology 1983, 23(3):345-357. 

33. Tomiyama K: Further observation on the time requirement for 
hypersensitive cell death of potatoes infected by Phytophthora infestans 
and its relation to metabolic activity. Phytopathologische Zeitschrift 1967, 
58(4):367-378. 

34. Kamoun S, Huitema E, VIeeshouwers V: Resistance to oomycetes: a general 
role for the hypersensitive response? Trends Plant Sci 1999, 4(5):1 96-200. 

35. Chen Y, Halterman DA: Phenotypic characterization of potato late blight 
resistance mediated by the broad-spectrum resistance gene RB. 
Phytopathology 201 1, 101 (2):263-270. 

36. Jung HW, Hwang BK: The leucine-rich repeat (LRR) protein, CaLRRI, 
interacts with the hypersensitive induced reaction (HIR) protein, CaHIRI, 
and suppresses cell death induced by the CaHIRI protein. Mol Plant 
Pathol 2007, 8(4):503-514. 

37. Jung HW, Lim CW, Lee SC, Choi HW, Hwang CH, Hwang BK: Distinct roles 
of the pepper hypersensitive induced reaction protein gene Cai-ilRl in 
disease and osmotic stress, as determined by comparative transcriptome 
and proteome analyses. Planta 2008, 227(2):409-425. 

38. Qi YP, Tsuda K, Nguyen LV, Wang X, Lin JS, Murphy AS, Glazebrook J, 
Thordal-Christensen H, Katagiri F: Physical association of arabidopsis 
hypersensitive induced reaction proteins (HIRs) with the immune 
receptor RPS2. J Biol Chem 201 1, 286(36):31 297-31 307. 

39. Yoshioka H, Sugie K, Park HJ, Maeda H, Tsuda N, Kawakita K, Doke N: 
Induction of plant gp91 phox homolog by fungal cell wall, arachidonic 
acid, and salicylic acid in potato. Mol Plant Microbe Interact 2001, 
14(6):725-736. 

40. Leon-Reyes A, Spoel SH, De Lange ES, Abe H, Kobayashi M, Tsuda S, 
Millenaar FF, Welschen RAM, Ritsema T, Pieterse CMJ: Ethylene modulates 
the role of nonexpressor of pathogenesis related Genesi in cross talk 
between salicylate and jasmonate signaling. Plant Physiol 2009, 
149(4):1 797-1 809. 

41. Niki T, Mitsuhara I, Seo S, Ohtsubo N, Ohashi Y: Antagonistic effect of 
salicylic acid and jasmonic acid on the expression of pathogenesis- 
related (PR) protein genes in wounded mature tobacco leaves. Plant Cell 
Physiol m8, 39(5):500-507. 

42. Nunez-Pastrana R, Arcos-Ortega GF, Souza-Perera RA, Sanchez-Borges CA, 
Nakazawa-Ueji YE, Garcia-Villalobos FJ, Guzman-Antonio AA, Zuniga-Aguilar 
JJ: Ethylene, but not salicylic acid or methyl jasmonate, induces a 
resistance response against Phytophthora capsici in habanero pepper. 
Eur J Plant Pathol 201 1 , 1 31 (4):669-683. 

43. Robert-Seilaniantz A, Grant M, Jones JDG: Hormone crosstalk in plant 
disease and defense: more than just jasmonate-salicylate antagonism. In 
Annual review of phytopathology Vol 49. Edited by VanAlfen NK, Bruening G, 
Leach JE. Palo Alto: Annual Reviews; 2011:317-343. 

44. Tumwine J, Frinking HD, Jeger MJ: Isolation techniques and cultural media 
for Phytophthora infestans from tomatoes. Mycologist 2000, 14(3):137-139. 

45. Cox M, Peterson D, Biggs P: SolexaQA: at-a-glance quality assessment of 
lllumina second-generation sequencing data. BMC Bioinforma 2010, 11(1):485. 

46. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, 
Salzberg SL, Rinn JL, Pachter L: Differential gene and transcript expression 



analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc 
2012, 7(3):562-578. 

47. Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient 
alignment of short DNA sequences to the human genome. Genome Biol 
2009, 10(3):R25. 

48. Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions 
with RNA-Seq. Bioinformatics 2009, 25(9):1 105-1 111. 

49. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, 
Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by 
RNA-Seq reveals unannotated transcripts and isoform switching during 
cell differentiation. Nat Biotechnol 2010, 28(5):51 1-U1 74. 

50. Benjamini Y, Hochberg Y: Controlling false discovery rate - a practical and 
powerful approach to multiple testing. J R Stat Soc Ser B-Methodol 1995, 
57(1):289-300. 

51. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display 
of genome-wide expression patterns. Proc Natl Acad Sci USA 1 998, 
95(25):1 4863-1 4868. 

52. Saldanha AJ: Java treeview-extensible visualization of microarray data. 
Bioinformatics 2004, 20(1 7):3246-3248. 

53. R Development Core Team: R: a language and environment for statistical 
computing. Vienna, Austria: R Foundation for Statistical Computing; 2012. 
URL http://www.R-project.org/. ISBN 3-900051-07-0. 

54. Thimm 0, Biasing 0, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller 
LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics 
data sets onto diagrams of metabolic pathways and other biological 
processes. Plant J 2004, 37(6):91 4-939. 

55. Millett BP, Bradeen JM: Development of allele-specific PCR and RT-PCR 
assays for clustered resistance genes using a potato late blight 
resistance transgene as a model. Theor AppI Genet 2007, 1 14(3):501-513. 

56. Nicot N, Hausman JF, Hoffmann L, Evers D: Housekeeping gene selection 
for real-time RT-PCR normalization in potato during biotic and abiotic 
stress. J Exp Bot 2005, 56(421 ):2907-291 4. 



doi:1 0.1 1 86/1 471 -21 64-1 4-340 

Cite this article as: Gao et al.: Insights into organ-specific pathogen 
defense responses in plants: RNA-seq analysis of potato tuber- 
Phytophthora infestans interactions. BMC Genomics 2013 14:340. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



